3.2.21 \(\int (1+c x)^2 (a+b \text {arctanh}(c x))^3 \, dx\) [121]

3.2.21.1 Optimal result
3.2.21.2 Mathematica [B] (verified)
3.2.21.3 Rubi [A] (verified)
3.2.21.4 Maple [C] (warning: unable to verify)
3.2.21.5 Fricas [F]
3.2.21.6 Sympy [F]
3.2.21.7 Maxima [F]
3.2.21.8 Giac [F]
3.2.21.9 Mupad [F(-1)]

3.2.21.1 Optimal result

Integrand size = 18, antiderivative size = 240 \[ \int (1+c x)^2 (a+b \text {arctanh}(c x))^3 \, dx=a b^2 x+b^3 x \text {arctanh}(c x)+\frac {5 b (a+b \text {arctanh}(c x))^2}{2 c}+3 b x (a+b \text {arctanh}(c x))^2+\frac {1}{2} b c x^2 (a+b \text {arctanh}(c x))^2+\frac {(1+c x)^3 (a+b \text {arctanh}(c x))^3}{3 c}-\frac {6 b^2 (a+b \text {arctanh}(c x)) \log \left (\frac {2}{1-c x}\right )}{c}-\frac {4 b (a+b \text {arctanh}(c x))^2 \log \left (\frac {2}{1-c x}\right )}{c}+\frac {b^3 \log \left (1-c^2 x^2\right )}{2 c}-\frac {3 b^3 \operatorname {PolyLog}\left (2,1-\frac {2}{1-c x}\right )}{c}-\frac {4 b^2 (a+b \text {arctanh}(c x)) \operatorname {PolyLog}\left (2,1-\frac {2}{1-c x}\right )}{c}+\frac {2 b^3 \operatorname {PolyLog}\left (3,1-\frac {2}{1-c x}\right )}{c} \]

output
a*b^2*x+b^3*x*arctanh(c*x)+5/2*b*(a+b*arctanh(c*x))^2/c+3*b*x*(a+b*arctanh 
(c*x))^2+1/2*b*c*x^2*(a+b*arctanh(c*x))^2+1/3*(c*x+1)^3*(a+b*arctanh(c*x)) 
^3/c-6*b^2*(a+b*arctanh(c*x))*ln(2/(-c*x+1))/c-4*b*(a+b*arctanh(c*x))^2*ln 
(2/(-c*x+1))/c+1/2*b^3*ln(-c^2*x^2+1)/c-3*b^3*polylog(2,1-2/(-c*x+1))/c-4* 
b^2*(a+b*arctanh(c*x))*polylog(2,1-2/(-c*x+1))/c+2*b^3*polylog(3,1-2/(-c*x 
+1))/c
 
3.2.21.2 Mathematica [B] (verified)

Leaf count is larger than twice the leaf count of optimal. \(488\) vs. \(2(240)=480\).

Time = 1.64 (sec) , antiderivative size = 488, normalized size of antiderivative = 2.03 \[ \int (1+c x)^2 (a+b \text {arctanh}(c x))^3 \, dx=\frac {6 a^3 c x+18 a^2 b c x+6 a b^2 c x+6 a^3 c^2 x^2+3 a^2 b c^2 x^2+2 a^3 c^3 x^3-6 a b^2 \text {arctanh}(c x)+18 a^2 b c x \text {arctanh}(c x)+36 a b^2 c x \text {arctanh}(c x)+6 b^3 c x \text {arctanh}(c x)+18 a^2 b c^2 x^2 \text {arctanh}(c x)+6 a b^2 c^2 x^2 \text {arctanh}(c x)+6 a^2 b c^3 x^3 \text {arctanh}(c x)-42 a b^2 \text {arctanh}(c x)^2-21 b^3 \text {arctanh}(c x)^2+18 a b^2 c x \text {arctanh}(c x)^2+18 b^3 c x \text {arctanh}(c x)^2+18 a b^2 c^2 x^2 \text {arctanh}(c x)^2+3 b^3 c^2 x^2 \text {arctanh}(c x)^2+6 a b^2 c^3 x^3 \text {arctanh}(c x)^2-14 b^3 \text {arctanh}(c x)^3+6 b^3 c x \text {arctanh}(c x)^3+6 b^3 c^2 x^2 \text {arctanh}(c x)^3+2 b^3 c^3 x^3 \text {arctanh}(c x)^3-48 a b^2 \text {arctanh}(c x) \log \left (1+e^{-2 \text {arctanh}(c x)}\right )-36 b^3 \text {arctanh}(c x) \log \left (1+e^{-2 \text {arctanh}(c x)}\right )-24 b^3 \text {arctanh}(c x)^2 \log \left (1+e^{-2 \text {arctanh}(c x)}\right )+21 a^2 b \log (1-c x)+3 a^2 b \log (1+c x)+18 a b^2 \log \left (1-c^2 x^2\right )+3 b^3 \log \left (1-c^2 x^2\right )+6 b^2 (4 a+3 b+4 b \text {arctanh}(c x)) \operatorname {PolyLog}\left (2,-e^{-2 \text {arctanh}(c x)}\right )+12 b^3 \operatorname {PolyLog}\left (3,-e^{-2 \text {arctanh}(c x)}\right )}{6 c} \]

input
Integrate[(1 + c*x)^2*(a + b*ArcTanh[c*x])^3,x]
 
output
(6*a^3*c*x + 18*a^2*b*c*x + 6*a*b^2*c*x + 6*a^3*c^2*x^2 + 3*a^2*b*c^2*x^2 
+ 2*a^3*c^3*x^3 - 6*a*b^2*ArcTanh[c*x] + 18*a^2*b*c*x*ArcTanh[c*x] + 36*a* 
b^2*c*x*ArcTanh[c*x] + 6*b^3*c*x*ArcTanh[c*x] + 18*a^2*b*c^2*x^2*ArcTanh[c 
*x] + 6*a*b^2*c^2*x^2*ArcTanh[c*x] + 6*a^2*b*c^3*x^3*ArcTanh[c*x] - 42*a*b 
^2*ArcTanh[c*x]^2 - 21*b^3*ArcTanh[c*x]^2 + 18*a*b^2*c*x*ArcTanh[c*x]^2 + 
18*b^3*c*x*ArcTanh[c*x]^2 + 18*a*b^2*c^2*x^2*ArcTanh[c*x]^2 + 3*b^3*c^2*x^ 
2*ArcTanh[c*x]^2 + 6*a*b^2*c^3*x^3*ArcTanh[c*x]^2 - 14*b^3*ArcTanh[c*x]^3 
+ 6*b^3*c*x*ArcTanh[c*x]^3 + 6*b^3*c^2*x^2*ArcTanh[c*x]^3 + 2*b^3*c^3*x^3* 
ArcTanh[c*x]^3 - 48*a*b^2*ArcTanh[c*x]*Log[1 + E^(-2*ArcTanh[c*x])] - 36*b 
^3*ArcTanh[c*x]*Log[1 + E^(-2*ArcTanh[c*x])] - 24*b^3*ArcTanh[c*x]^2*Log[1 
 + E^(-2*ArcTanh[c*x])] + 21*a^2*b*Log[1 - c*x] + 3*a^2*b*Log[1 + c*x] + 1 
8*a*b^2*Log[1 - c^2*x^2] + 3*b^3*Log[1 - c^2*x^2] + 6*b^2*(4*a + 3*b + 4*b 
*ArcTanh[c*x])*PolyLog[2, -E^(-2*ArcTanh[c*x])] + 12*b^3*PolyLog[3, -E^(-2 
*ArcTanh[c*x])])/(6*c)
 
3.2.21.3 Rubi [A] (verified)

Time = 0.69 (sec) , antiderivative size = 236, normalized size of antiderivative = 0.98, number of steps used = 2, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.111, Rules used = {6480, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int (c x+1)^2 (a+b \text {arctanh}(c x))^3 \, dx\)

\(\Big \downarrow \) 6480

\(\displaystyle \frac {(c x+1)^3 (a+b \text {arctanh}(c x))^3}{3 c}-b \int \left (-c x (a+b \text {arctanh}(c x))^2+\frac {4 (c x+1) (a+b \text {arctanh}(c x))^2}{1-c^2 x^2}-3 (a+b \text {arctanh}(c x))^2\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {(c x+1)^3 (a+b \text {arctanh}(c x))^3}{3 c}-b \left (\frac {4 b \operatorname {PolyLog}\left (2,1-\frac {2}{1-c x}\right ) (a+b \text {arctanh}(c x))}{c}-\frac {1}{2} c x^2 (a+b \text {arctanh}(c x))^2-3 x (a+b \text {arctanh}(c x))^2-\frac {5 (a+b \text {arctanh}(c x))^2}{2 c}+\frac {6 b \log \left (\frac {2}{1-c x}\right ) (a+b \text {arctanh}(c x))}{c}+\frac {4 \log \left (\frac {2}{1-c x}\right ) (a+b \text {arctanh}(c x))^2}{c}-a b x+b^2 (-x) \text {arctanh}(c x)-\frac {b^2 \log \left (1-c^2 x^2\right )}{2 c}+\frac {3 b^2 \operatorname {PolyLog}\left (2,1-\frac {2}{1-c x}\right )}{c}-\frac {2 b^2 \operatorname {PolyLog}\left (3,1-\frac {2}{1-c x}\right )}{c}\right )\)

input
Int[(1 + c*x)^2*(a + b*ArcTanh[c*x])^3,x]
 
output
((1 + c*x)^3*(a + b*ArcTanh[c*x])^3)/(3*c) - b*(-(a*b*x) - b^2*x*ArcTanh[c 
*x] - (5*(a + b*ArcTanh[c*x])^2)/(2*c) - 3*x*(a + b*ArcTanh[c*x])^2 - (c*x 
^2*(a + b*ArcTanh[c*x])^2)/2 + (6*b*(a + b*ArcTanh[c*x])*Log[2/(1 - c*x)]) 
/c + (4*(a + b*ArcTanh[c*x])^2*Log[2/(1 - c*x)])/c - (b^2*Log[1 - c^2*x^2] 
)/(2*c) + (3*b^2*PolyLog[2, 1 - 2/(1 - c*x)])/c + (4*b*(a + b*ArcTanh[c*x] 
)*PolyLog[2, 1 - 2/(1 - c*x)])/c - (2*b^2*PolyLog[3, 1 - 2/(1 - c*x)])/c)
 

3.2.21.3.1 Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 6480
Int[((a_.) + ArcTanh[(c_.)*(x_)]*(b_.))^(p_)*((d_) + (e_.)*(x_))^(q_.), x_S 
ymbol] :> Simp[(d + e*x)^(q + 1)*((a + b*ArcTanh[c*x])^p/(e*(q + 1))), x] - 
 Simp[b*c*(p/(e*(q + 1)))   Int[ExpandIntegrand[(a + b*ArcTanh[c*x])^(p - 1 
), (d + e*x)^(q + 1)/(1 - c^2*x^2), x], x], x] /; FreeQ[{a, b, c, d, e}, x] 
 && IGtQ[p, 1] && IntegerQ[q] && NeQ[q, -1]
 
3.2.21.4 Maple [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 1.87 (sec) , antiderivative size = 610, normalized size of antiderivative = 2.54

method result size
derivativedivides \(\frac {\frac {\left (c x +1\right )^{3} a^{3}}{3}+b^{3} \left (\frac {\operatorname {arctanh}\left (c x \right )^{3} c^{3} x^{3}}{3}+\operatorname {arctanh}\left (c x \right )^{3} c^{2} x^{2}+\operatorname {arctanh}\left (c x \right )^{3} c x +\frac {\operatorname {arctanh}\left (c x \right )^{3}}{3}+\frac {c^{2} x^{2} \operatorname {arctanh}\left (c x \right )^{2}}{2}+3 c x \operatorname {arctanh}\left (c x \right )^{2}+4 \operatorname {arctanh}\left (c x \right )^{2} \ln \left (c x -1\right )-4 \,\operatorname {arctanh}\left (c x \right ) \operatorname {polylog}\left (2, -\frac {\left (c x +1\right )^{2}}{-c^{2} x^{2}+1}\right )+2 \operatorname {polylog}\left (3, -\frac {\left (c x +1\right )^{2}}{-c^{2} x^{2}+1}\right )-\ln \left (1+\frac {\left (c x +1\right )^{2}}{-c^{2} x^{2}+1}\right )+\frac {5 \operatorname {arctanh}\left (c x \right )^{2}}{2}-6 \operatorname {dilog}\left (1+\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )-6 \operatorname {dilog}\left (1-\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )-6 \,\operatorname {arctanh}\left (c x \right ) \ln \left (1+\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )-6 \,\operatorname {arctanh}\left (c x \right ) \ln \left (1-\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )+\left (c x +1\right ) \operatorname {arctanh}\left (c x \right )+4 i \pi {\operatorname {csgn}\left (\frac {i}{1-\frac {\left (c x +1\right )^{2}}{c^{2} x^{2}-1}}\right )}^{2} \operatorname {arctanh}\left (c x \right )^{2}-4 \operatorname {arctanh}\left (c x \right )^{2} \ln \left (2\right )-4 i \pi {\operatorname {csgn}\left (\frac {i}{1-\frac {\left (c x +1\right )^{2}}{c^{2} x^{2}-1}}\right )}^{3} \operatorname {arctanh}\left (c x \right )^{2}-4 i \operatorname {arctanh}\left (c x \right )^{2} \pi \right )+3 a \,b^{2} \left (\frac {\operatorname {arctanh}\left (c x \right )^{2} c^{3} x^{3}}{3}+c^{2} x^{2} \operatorname {arctanh}\left (c x \right )^{2}+c x \operatorname {arctanh}\left (c x \right )^{2}+\frac {\operatorname {arctanh}\left (c x \right )^{2}}{3}+\frac {c^{2} x^{2} \operatorname {arctanh}\left (c x \right )}{3}+2 c x \,\operatorname {arctanh}\left (c x \right )+\frac {8 \,\operatorname {arctanh}\left (c x \right ) \ln \left (c x -1\right )}{3}+\frac {c x}{3}-\frac {1}{3}+\frac {5 \ln \left (c x +1\right )}{6}+\frac {7 \ln \left (c x -1\right )}{6}-\frac {4 \operatorname {dilog}\left (\frac {c x}{2}+\frac {1}{2}\right )}{3}-\frac {4 \ln \left (c x -1\right ) \ln \left (\frac {c x}{2}+\frac {1}{2}\right )}{3}+\frac {2 \ln \left (c x -1\right )^{2}}{3}\right )+3 a^{2} b \left (\frac {c^{3} x^{3} \operatorname {arctanh}\left (c x \right )}{3}+c^{2} x^{2} \operatorname {arctanh}\left (c x \right )+c x \,\operatorname {arctanh}\left (c x \right )+\frac {\operatorname {arctanh}\left (c x \right )}{3}+\frac {c^{2} x^{2}}{6}+c x +\frac {4 \ln \left (c x -1\right )}{3}\right )}{c}\) \(610\)
default \(\frac {\frac {\left (c x +1\right )^{3} a^{3}}{3}+b^{3} \left (\frac {\operatorname {arctanh}\left (c x \right )^{3} c^{3} x^{3}}{3}+\operatorname {arctanh}\left (c x \right )^{3} c^{2} x^{2}+\operatorname {arctanh}\left (c x \right )^{3} c x +\frac {\operatorname {arctanh}\left (c x \right )^{3}}{3}+\frac {c^{2} x^{2} \operatorname {arctanh}\left (c x \right )^{2}}{2}+3 c x \operatorname {arctanh}\left (c x \right )^{2}+4 \operatorname {arctanh}\left (c x \right )^{2} \ln \left (c x -1\right )-4 \,\operatorname {arctanh}\left (c x \right ) \operatorname {polylog}\left (2, -\frac {\left (c x +1\right )^{2}}{-c^{2} x^{2}+1}\right )+2 \operatorname {polylog}\left (3, -\frac {\left (c x +1\right )^{2}}{-c^{2} x^{2}+1}\right )-\ln \left (1+\frac {\left (c x +1\right )^{2}}{-c^{2} x^{2}+1}\right )+\frac {5 \operatorname {arctanh}\left (c x \right )^{2}}{2}-6 \operatorname {dilog}\left (1+\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )-6 \operatorname {dilog}\left (1-\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )-6 \,\operatorname {arctanh}\left (c x \right ) \ln \left (1+\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )-6 \,\operatorname {arctanh}\left (c x \right ) \ln \left (1-\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )+\left (c x +1\right ) \operatorname {arctanh}\left (c x \right )+4 i \pi {\operatorname {csgn}\left (\frac {i}{1-\frac {\left (c x +1\right )^{2}}{c^{2} x^{2}-1}}\right )}^{2} \operatorname {arctanh}\left (c x \right )^{2}-4 \operatorname {arctanh}\left (c x \right )^{2} \ln \left (2\right )-4 i \pi {\operatorname {csgn}\left (\frac {i}{1-\frac {\left (c x +1\right )^{2}}{c^{2} x^{2}-1}}\right )}^{3} \operatorname {arctanh}\left (c x \right )^{2}-4 i \operatorname {arctanh}\left (c x \right )^{2} \pi \right )+3 a \,b^{2} \left (\frac {\operatorname {arctanh}\left (c x \right )^{2} c^{3} x^{3}}{3}+c^{2} x^{2} \operatorname {arctanh}\left (c x \right )^{2}+c x \operatorname {arctanh}\left (c x \right )^{2}+\frac {\operatorname {arctanh}\left (c x \right )^{2}}{3}+\frac {c^{2} x^{2} \operatorname {arctanh}\left (c x \right )}{3}+2 c x \,\operatorname {arctanh}\left (c x \right )+\frac {8 \,\operatorname {arctanh}\left (c x \right ) \ln \left (c x -1\right )}{3}+\frac {c x}{3}-\frac {1}{3}+\frac {5 \ln \left (c x +1\right )}{6}+\frac {7 \ln \left (c x -1\right )}{6}-\frac {4 \operatorname {dilog}\left (\frac {c x}{2}+\frac {1}{2}\right )}{3}-\frac {4 \ln \left (c x -1\right ) \ln \left (\frac {c x}{2}+\frac {1}{2}\right )}{3}+\frac {2 \ln \left (c x -1\right )^{2}}{3}\right )+3 a^{2} b \left (\frac {c^{3} x^{3} \operatorname {arctanh}\left (c x \right )}{3}+c^{2} x^{2} \operatorname {arctanh}\left (c x \right )+c x \,\operatorname {arctanh}\left (c x \right )+\frac {\operatorname {arctanh}\left (c x \right )}{3}+\frac {c^{2} x^{2}}{6}+c x +\frac {4 \ln \left (c x -1\right )}{3}\right )}{c}\) \(610\)
parts \(\frac {a^{3} \left (c x +1\right )^{3}}{3 c}+\frac {b^{3} \left (\frac {\operatorname {arctanh}\left (c x \right )^{3} c^{3} x^{3}}{3}+\operatorname {arctanh}\left (c x \right )^{3} c^{2} x^{2}+\operatorname {arctanh}\left (c x \right )^{3} c x +\frac {\operatorname {arctanh}\left (c x \right )^{3}}{3}+\frac {c^{2} x^{2} \operatorname {arctanh}\left (c x \right )^{2}}{2}+3 c x \operatorname {arctanh}\left (c x \right )^{2}+4 \operatorname {arctanh}\left (c x \right )^{2} \ln \left (c x -1\right )-4 \,\operatorname {arctanh}\left (c x \right ) \operatorname {polylog}\left (2, -\frac {\left (c x +1\right )^{2}}{-c^{2} x^{2}+1}\right )+2 \operatorname {polylog}\left (3, -\frac {\left (c x +1\right )^{2}}{-c^{2} x^{2}+1}\right )-\ln \left (1+\frac {\left (c x +1\right )^{2}}{-c^{2} x^{2}+1}\right )+\frac {5 \operatorname {arctanh}\left (c x \right )^{2}}{2}-6 \operatorname {dilog}\left (1+\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )-6 \operatorname {dilog}\left (1-\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )-6 \,\operatorname {arctanh}\left (c x \right ) \ln \left (1+\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )-6 \,\operatorname {arctanh}\left (c x \right ) \ln \left (1-\frac {i \left (c x +1\right )}{\sqrt {-c^{2} x^{2}+1}}\right )+\left (c x +1\right ) \operatorname {arctanh}\left (c x \right )+4 i \pi {\operatorname {csgn}\left (\frac {i}{1-\frac {\left (c x +1\right )^{2}}{c^{2} x^{2}-1}}\right )}^{2} \operatorname {arctanh}\left (c x \right )^{2}-4 \operatorname {arctanh}\left (c x \right )^{2} \ln \left (2\right )-4 i \pi {\operatorname {csgn}\left (\frac {i}{1-\frac {\left (c x +1\right )^{2}}{c^{2} x^{2}-1}}\right )}^{3} \operatorname {arctanh}\left (c x \right )^{2}-4 i \operatorname {arctanh}\left (c x \right )^{2} \pi \right )}{c}+\frac {3 a \,b^{2} \left (\frac {\operatorname {arctanh}\left (c x \right )^{2} c^{3} x^{3}}{3}+c^{2} x^{2} \operatorname {arctanh}\left (c x \right )^{2}+c x \operatorname {arctanh}\left (c x \right )^{2}+\frac {\operatorname {arctanh}\left (c x \right )^{2}}{3}+\frac {c^{2} x^{2} \operatorname {arctanh}\left (c x \right )}{3}+2 c x \,\operatorname {arctanh}\left (c x \right )+\frac {8 \,\operatorname {arctanh}\left (c x \right ) \ln \left (c x -1\right )}{3}+\frac {c x}{3}-\frac {1}{3}+\frac {5 \ln \left (c x +1\right )}{6}+\frac {7 \ln \left (c x -1\right )}{6}-\frac {4 \operatorname {dilog}\left (\frac {c x}{2}+\frac {1}{2}\right )}{3}-\frac {4 \ln \left (c x -1\right ) \ln \left (\frac {c x}{2}+\frac {1}{2}\right )}{3}+\frac {2 \ln \left (c x -1\right )^{2}}{3}\right )}{c}+\frac {3 a^{2} b \left (\frac {c^{3} x^{3} \operatorname {arctanh}\left (c x \right )}{3}+c^{2} x^{2} \operatorname {arctanh}\left (c x \right )+c x \,\operatorname {arctanh}\left (c x \right )+\frac {\operatorname {arctanh}\left (c x \right )}{3}+\frac {c^{2} x^{2}}{6}+c x +\frac {4 \ln \left (c x -1\right )}{3}\right )}{c}\) \(618\)

input
int((c*x+1)^2*(a+b*arctanh(c*x))^3,x,method=_RETURNVERBOSE)
 
output
1/c*(1/3*(c*x+1)^3*a^3+b^3*(1/3*arctanh(c*x)^3*c^3*x^3+arctanh(c*x)^3*c^2* 
x^2+arctanh(c*x)^3*c*x+1/3*arctanh(c*x)^3+1/2*c^2*x^2*arctanh(c*x)^2+3*c*x 
*arctanh(c*x)^2+4*arctanh(c*x)^2*ln(c*x-1)-4*arctanh(c*x)*polylog(2,-(c*x+ 
1)^2/(-c^2*x^2+1))+2*polylog(3,-(c*x+1)^2/(-c^2*x^2+1))-ln(1+(c*x+1)^2/(-c 
^2*x^2+1))+5/2*arctanh(c*x)^2-6*dilog(1+I*(c*x+1)/(-c^2*x^2+1)^(1/2))-6*di 
log(1-I*(c*x+1)/(-c^2*x^2+1)^(1/2))-6*arctanh(c*x)*ln(1+I*(c*x+1)/(-c^2*x^ 
2+1)^(1/2))-6*arctanh(c*x)*ln(1-I*(c*x+1)/(-c^2*x^2+1)^(1/2))+(c*x+1)*arct 
anh(c*x)+4*I*Pi*csgn(I/(1-(c*x+1)^2/(c^2*x^2-1)))^2*arctanh(c*x)^2-4*arcta 
nh(c*x)^2*ln(2)-4*I*Pi*csgn(I/(1-(c*x+1)^2/(c^2*x^2-1)))^3*arctanh(c*x)^2- 
4*I*arctanh(c*x)^2*Pi)+3*a*b^2*(1/3*arctanh(c*x)^2*c^3*x^3+c^2*x^2*arctanh 
(c*x)^2+c*x*arctanh(c*x)^2+1/3*arctanh(c*x)^2+1/3*c^2*x^2*arctanh(c*x)+2*c 
*x*arctanh(c*x)+8/3*arctanh(c*x)*ln(c*x-1)+1/3*c*x-1/3+5/6*ln(c*x+1)+7/6*l 
n(c*x-1)-4/3*dilog(1/2*c*x+1/2)-4/3*ln(c*x-1)*ln(1/2*c*x+1/2)+2/3*ln(c*x-1 
)^2)+3*a^2*b*(1/3*c^3*x^3*arctanh(c*x)+c^2*x^2*arctanh(c*x)+c*x*arctanh(c* 
x)+1/3*arctanh(c*x)+1/6*c^2*x^2+c*x+4/3*ln(c*x-1)))
 
3.2.21.5 Fricas [F]

\[ \int (1+c x)^2 (a+b \text {arctanh}(c x))^3 \, dx=\int { {\left (c x + 1\right )}^{2} {\left (b \operatorname {artanh}\left (c x\right ) + a\right )}^{3} \,d x } \]

input
integrate((c*x+1)^2*(a+b*arctanh(c*x))^3,x, algorithm="fricas")
 
output
integral(a^3*c^2*x^2 + 2*a^3*c*x + (b^3*c^2*x^2 + 2*b^3*c*x + b^3)*arctanh 
(c*x)^3 + a^3 + 3*(a*b^2*c^2*x^2 + 2*a*b^2*c*x + a*b^2)*arctanh(c*x)^2 + 3 
*(a^2*b*c^2*x^2 + 2*a^2*b*c*x + a^2*b)*arctanh(c*x), x)
 
3.2.21.6 Sympy [F]

\[ \int (1+c x)^2 (a+b \text {arctanh}(c x))^3 \, dx=\int \left (a + b \operatorname {atanh}{\left (c x \right )}\right )^{3} \left (c x + 1\right )^{2}\, dx \]

input
integrate((c*x+1)**2*(a+b*atanh(c*x))**3,x)
 
output
Integral((a + b*atanh(c*x))**3*(c*x + 1)**2, x)
 
3.2.21.7 Maxima [F]

\[ \int (1+c x)^2 (a+b \text {arctanh}(c x))^3 \, dx=\int { {\left (c x + 1\right )}^{2} {\left (b \operatorname {artanh}\left (c x\right ) + a\right )}^{3} \,d x } \]

input
integrate((c*x+1)^2*(a+b*arctanh(c*x))^3,x, algorithm="maxima")
 
output
1/3*a^3*c^2*x^3 + 1/2*(2*x^3*arctanh(c*x) + c*(x^2/c^2 + log(c^2*x^2 - 1)/ 
c^4))*a^2*b*c^2 + a^3*c*x^2 + 3/2*(2*x^2*arctanh(c*x) + c*(2*x/c^2 - log(c 
*x + 1)/c^3 + log(c*x - 1)/c^3))*a^2*b*c + a^3*x + 3/2*(2*c*x*arctanh(c*x) 
 + log(-c^2*x^2 + 1))*a^2*b/c - 1/24*((b^3*c^3*x^3 + 3*b^3*c^2*x^2 + 3*b^3 
*c*x - 7*b^3)*log(-c*x + 1)^3 - 3*(2*a*b^2*c^3*x^3 + (6*a*b^2*c^2 + b^3*c^ 
2)*x^2 + 6*(a*b^2*c + b^3*c)*x + (b^3*c^3*x^3 + 3*b^3*c^2*x^2 + 3*b^3*c*x 
+ b^3)*log(c*x + 1))*log(-c*x + 1)^2)/c - integrate(-1/8*((b^3*c^3*x^3 + b 
^3*c^2*x^2 - b^3*c*x - b^3)*log(c*x + 1)^3 + 6*(a*b^2*c^3*x^3 + a*b^2*c^2* 
x^2 - a*b^2*c*x - a*b^2)*log(c*x + 1)^2 - (4*a*b^2*c^3*x^3 + 2*(6*a*b^2*c^ 
2 + b^3*c^2)*x^2 + 3*(b^3*c^3*x^3 + b^3*c^2*x^2 - b^3*c*x - b^3)*log(c*x + 
 1)^2 + 12*(a*b^2*c + b^3*c)*x + 2*((6*a*b^2*c^3 + b^3*c^3)*x^3 - 6*a*b^2 
+ b^3 + 3*(2*a*b^2*c^2 + b^3*c^2)*x^2 - 3*(2*a*b^2*c - b^3*c)*x)*log(c*x + 
 1))*log(-c*x + 1))/(c*x - 1), x)
 
3.2.21.8 Giac [F]

\[ \int (1+c x)^2 (a+b \text {arctanh}(c x))^3 \, dx=\int { {\left (c x + 1\right )}^{2} {\left (b \operatorname {artanh}\left (c x\right ) + a\right )}^{3} \,d x } \]

input
integrate((c*x+1)^2*(a+b*arctanh(c*x))^3,x, algorithm="giac")
 
output
integrate((c*x + 1)^2*(b*arctanh(c*x) + a)^3, x)
 
3.2.21.9 Mupad [F(-1)]

Timed out. \[ \int (1+c x)^2 (a+b \text {arctanh}(c x))^3 \, dx=\int {\left (a+b\,\mathrm {atanh}\left (c\,x\right )\right )}^3\,{\left (c\,x+1\right )}^2 \,d x \]

input
int((a + b*atanh(c*x))^3*(c*x + 1)^2,x)
 
output
int((a + b*atanh(c*x))^3*(c*x + 1)^2, x)